Universal occurrence of mixed-synchronization in counter-rotating nonlinear coupled 

oscillators 
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By coupling counter-rotating coupled nonlinear oscillators, we observe a "mixed" synchronization 
between the different dynamical variables of the same system. The phenomenon of amplitude death 
is also observed. Results for coupled systems with co-rotating coupled oscillators are also presented 
for a detailed comparison. Results for Landau-Stuart and Rossler oscillators are presented. 
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Natural systems are rarely isolated, and thus studies 
of coupled dynamical systems which arise in a variety of 
contexts in the physical, biological, and social sciences, 
have broad relevance to many areas of research. 

Consider the coupling of two nonlinear oscillators, 
schematically shown in Fig. [T] with the symbols and 
'-' indicating the sense of rotation of the orbits in the 
phase space, namely clockwise or counter-clockwise. If 
directions of rotation of both oscillators are the same, 
as in Fig. [TJa), the system is co- rotating: such coupled 
nonlinear oscillators have been extensively studied from 
both theoretical and experimental points of view [l| . 

Important phenomena which have been observed for 
co-rotating oscillators include synchronization, hystere- 
hase locking, phase shifting, phase-flip, or riddling 
There are many forms of synchrony e. g. com- 
lete synchronization: system variables become identical 
phase synchronization: phase difference between the 
oscillators becomes bounded 0], and others e. g. lag syn- 
chronization and generalized synchronization etc. In 
all these type of synchronizations if one of the dynamical 
variable is synchronized then rest of the variables follow 
the same. 

Another important property of co-rotating oscillators 
is amplitude death (AD) [8( where amplitude of the oscil- 
lation ceases to zero at stable fixed point of the system. 
This AD can be achieved with various type of interac- 
tions between the oscillators e. g. mismatched oscillators 
@, delayed [l(|, conjugate [ll| and nonlinear [l2| etc. 

When the rotations of the individual oscillators are in 
opposite senses, as in Fig. Hlb), then the system is termed 
counter-rotating. To the best of our knowledge, counter- 
rotating systems have not explicitly been studied earlier, 
and in this Letter, we focus on such systems. We ob- 
serve that both AD as well as synchronization arise with 
this form of coupling. A novel "mixed" type of synchro- 
nization among different variables is seen to occur: some 
variables are synchronized in-phase while other variables 
can be out-of-phase. 




FIG. I: Schematic coupling schemes for two coupled oscilla- 
tors: (a) co-rotating and (b) counter-rotating. Arrows show 
the directions of rotation while — and + indicate the corre- 
sponding clockwise and anticlockwise rotations respectively. 
H denotes the coupling function, Eq. |T}. 



Consider a coupled system of N nonlinear oscillators 
specified by the equations 



X, = F,(X 4> u; i ) 



^^A^HfJ^X, 

1 3 = 1 



.,N(1) 



where X, is m-dimensional vector of dynamical variables 
and the Fj's specify the evolution equations of the dif- 
ferent oscillators with internal frequency parameter a;,. 
Here Ki is the number of connections to the i— th oscil- 
lators, namely its degree such that 1 < Ki < N, and 
e is the coupling strength. The connection topology is 
specified by the adjacency matrix A, with Aij — 1 if os- 
cillators i and j are connected to each other and Ay = 
otherwise. The coupling function H : W 71 — > M. m specifies 
the manner in which the oscillators i and j are coupled, 
with H(Xj, Xj) being a function of Xj and Xj. 

We first analyze two Landau-Stuart limit cycle oscilla- 
tors [6| coupled in the co- (Fig. [2(a)) as well as counter- 
rotating (Fig. Hlb)) schemes. The coupled equations are 



Zi = (1 



- \Zi\ 2 )Zi+e(Zj -Z t ) 



(2) 
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where i,j = 1,2, i ^ j and Zi is the complex variables 
with real part Xi and imaginary part yi. In Eq. ([1]) the 
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variables are X t = [x u y t ] T and H = [(Xj—Xi), {y 3 -yi)] T 
where T denotes the transpose. In the uncoupled case, 



0, each oscillator has a fixed point at 







i.e. e 

which is unstable with eigenvalues 1 ± iu)j . Irrespective 
of the values of u>j each oscillator has the same attractor 
with unit radius, and corresponding frequency ujj . Shown 
in Fig. H](a) and (b) are the transient trajectories for co- 
{uji = LU2 — 5) and counter-rotating (wi = 5,lu2 — —5) 
respectively for two uncoupled oscillators. These clearly 
show that, in these two cases, the trajectories are moving 
in the same or in the opposite directions respectively. In 
order to see the effects of coupling on these both schemes, 
results for co- and counter rotating coupled oscillators 
are shown in the left and the right panels respectively 
in Fig. [2] As the coupling is introduced, the dynami- 
cal behavior changes. The few largest Lyapunov expo- 
nents (see Fig. [UJc,d)) clearly show that in the case of 
co-rotating oscillators, the largest Lyapunov exponent re- 
mains zero for all coupling strengths which suggests the 
continuation of periodic motion. This is shown Fig. [2Je) 
at e = 0.5. However for counter-rotating oscillators, the 
maximal Lyapunov exponent (Fig. [2Id)) becomes neg- 
ative (this is an instance of amplitude death, marked 
"AD") where the fixed point Z*x = Z„ 2 = is stabi- 
lized. A simple eigenvalue (A) calculation, with real part 
a and imaginary part /3, gives the relation, 



a 



1 — e ± e and (3 



(3) 



where u>\ = uji = oj for co-rotation (left panel). For co- 
rotating coupled oscillators the real part of the largest 
eigenvalue, a = 1, and hence fixed point will be unstable 
for all values of oj and e. This rules out the possibility 
of AD in such scheme (for details see Ref. Q). However 
for uj\ ~ uj,lo2 = —lu for the case of counter-rotating 
oscillators (right panel) the eigenvalues will be 



a = 1 — e and (3 = \J uj 2 — e 2 
for e < oj 2 , while 

a = 1 — e ± \J e 2 — uj 2 and (3 = 0, 



(4) 



(5) 



for e 2 > uj 2 . These relations give the region 1 < e < 
(1 + lu 2 )/2 where the fixed point at the origin will be 
stable (a < 0) as shown in Fig. [2fd) by the arrow. The 
transient trajectory in AD region at e = 2 is shown in 
the inset figure. This clearly shows, which is a new re- 
sult, that AD is possible (via a Hopf bifurcation) in the 
counter-rotating identical oscillators while it is impossi- 
ble in the co-rotating case. It should be noted that AD is 
possible in co-rotating case only when there are mismatch 
in frequencies of the oscillators 9( . Here we show that in 
counter-rotating case everything is same as that of the co- 
rotating except directions of rotation. This novel result 
clearly distinguishes counter-rotating oscillators from the 
more extensively studied co-rotating case. 

Another important result pertains to the regime of syn- 
chronized motion. Figure[5Je) shows the attractors which 
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FIG. 2: Left and right panels correspond to the co- and 
counter-rotating respectively Landau-Stuart oscillators, Eq. 
(J2J. (a,b) The transient trajectories of uncoupled oscilla- 
tors. (c,d) A few largest Lyapunov exponents with coupling 
strength, e. Inset figure in (d) corresponds to the transient 
trajectories in AD regime at e = 2. (e,f) The asymptotic at- 
tractors in Xi — yi plane. (g,h) The signature of relative phase 
in x\ vs X2 and j/i vs yi at e = 0.5. Dashed lines in (h) show 
the relative phase-difference of (black) and 7r (red). 



are exactly same as uncoupled case (Fig. Ufa)) having 
complete synchronization. The relative phase between 
the oscillators is shown in Fig. [UJg) where trajectories 
lie on the synchronization manifold, x\ = X2,y\ = y2) 
this confirms the occurrence of complete synchronization. 
However in the counter-rotating case, Fig. [2jf ) the at- 
tractors are not identical: they do not remain as limit cy- 
cle circles (with unit radius) instead change shape as well 
as get reduced in amplitude. The oscillators are synchro- 
nized though, as can be seen in the relative phase shown 
in Fig. [Hff). The variables essentially in-phase while 
y— variables are out-of-phase (dashed lines correspond 
to the relative phase zero (black) and ir (red)). The si- 
multaneous presence of nearly in- as well as out-of-phase 
synchrony in different variables among the same oscilla- 
tors is termed here as "mixed" synchronization. 

It is important to note that the usual definition of 
phase difference between oscillators is not applicable here 
since there is continuous change of individual angle in 
opposite directions. Therefore the effect of counter- 
rotation has drastic difference due to presence of mixed- 
synchronization as compared to co-rotation. 

The generality of these results is demonstrated in a 
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FIG. 3: (a) The largest two Lyapunov exponents with in- 
ternal parameter, for anticlockwise (cui = 1 - black lines) 
and clockwise (uj2 = — 1 - red lines) for uncoupled (e = 0) 
Rossler oscillators, Eq. (|6]). The trajectories of individual 
oscillators for (b) periodic (fi = 2) and (c) chaotic (/i = 10) 
motions. 



system of coupled Rossler oscillators [13|, given by the 
equations of motion 



dxj (t) 

dt 
dyijt) 

dt 
dzj(t) 
dt 



-uJiVi - Zi + e(xj - Xi) 
uj % Xi + Q.2j/j 
0.2 + Zi(xi - fi) 



(6) 



where u>i (internal frequency) and /i are parameters, 
and i,j — 1,2, i ^ j. (With regard to Eq. (JTJ, 
Xi = [x il y i ,z i ] T and H = [(x - Xj),0,0] T .) 

Changing the parameter u>i merely changes the sense 
of rotation; as can be seen in Fig. [3ja) the two largest 
Lyapunov exponents as a function of fi are essentially 
identical for u>\ = 1 (black curves) and ll>2 = — 1 (red 
curves). Corresponding typical trajectories for periodic 
and chaotic motions are shown in Fig. ISJb) and (c) at 
fi = 2 and 10 respectively. Note that there is change of 
position of the attractors in the phase space. These fixed 
points (x*i = (fi± ^ n 2 - 0.16/o; l 2 )/2, y H = 5ujiX ti ,z„ = 
5w^a;*i) depend on the sign of uii (unlike the case of Eq. 
(|2"|) where the origin fixed point is independent of fre- 
quency) . 

The change in dynamical behavior arising from the 
coupling is shown in Fig. |4] for co- (left panel) and 
counter-rotating (right panel) coupled chaotic (fi = 10) 
Rossler oscillators, Eq. JBJ. For small coupling, say 
e = 0.005, the motion is unsynchronized which can be 
seen from the relative phase between the trajectories in 
Fig. |U(c) in x\ vs X2 (black) and y\ vs 7/2 (red). As the 
coupling strength is increased beyond e ~ 0.02 phase- 
synchronization sets in (the fourth largest Lyapunov ex- 
ponent (in green) becomes negative [3]), and after this all 
variables x, y, z remain in-phase as shown in Fig. IHe). 
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FIG. 4: Left and right panels correspond to the co- and 
counter-rotating chaotic Rossler oscillators respectively, Eq. 
§6§. (a,b) A few largest Lyapunov exponents with cou- 
pling strength, e at fi = 10. The signature of (c,d) 
non-synchronization, (e) phase-synchronization, (f) mixed- 
synchronization, (g) complete-synchronization and (h) mixed- 
synchronization. The coupling strengths for (c,d), (e,f) and 
(g,h) are at e = 0.005, 0.05 and 0.2 respectively. 



When the coupling is increased beyond e ~ 0.12 com- 
plete synchronization occurs as seen in Fig. @] (g) 
where x x = x 2 ,yi = yi- 

For counter-rotating oscillators (ujx = l,o>2 = — 1) 
(right panel) Lyapunov spectrum (Fig. 0Jb)) remains 
the same as for co-rotating oscillators (Fig. Ufa)) show- 
ing that the invariant measure does not change. For small 
coupling the oscillators are not in the phase synchronized 
state Fig. EJd), as analogous to the co-rotating case Fig. 
SJc). However, as coupling is increased further, the phase 
synchronization sets in, (Fig. H)(f)), but the relative phase 
between the three variables is not identical. 

Shown in Fig. [4] (f) is the relative phase in x\ — X2 
and 2/1 — 2/2 planes. Clearly x\ and xi are in-phase 
(black) synchronized while y\ and 2/2 are out-of-phase 
(red). This is clearer for higher coupling, e = 0.2 (Fig. 
|4jh)) where there is complete synchronization in X\ and 
X2 (zero relative phase) while y\ and 2/2 remain out-of- 
phase with phase difference it. We find z\ and Z2 in 
complete synchrony (results not shown here (Tij). For 
counter-rotating coupled oscillators it appears that if one 
variable is in-phase synchrony, then at least one other 
variable must be out-of-phase, i.e. if motion is on the 
manifold {x\ = X2,y\ — 2/2,^1 = Z2) for co-rotating 
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FIG. 5: Left panel shows the network of different topolo- 
gies for TV = 10 coupled oscillators: (a) periodic chain (b) 
globally coupled and (c) small-world network-see the text 
for details. The largest three Lyapunov exponents with cou- 
pling strength for different topologies are plotted in the middle 
and right panel corresponding to the Landau-Stuart and the 
Rossler coupled oscillators respectively. Inset figures show the 
relative phase of xa vs xb (black) and yA vs ys (red) between 
the marked oscillators A and B. The coupling strengths for 
the inset figures of middle panels (d,e,f) are at e = 1.5 while 
the corresponding values for (g), (h) and (i) are at e = 2,0.5 
and 1 respectively. 



case (Fig. S[g)) then the corresponding manifold for the 
counter-rotating case is [x\ — £2,2/1 = — 2/2>^i = z 2), 
Fig. Hh). 

Mixed-synchronization and AD also occur in network 
of oscillators where the rotation sense is chosen randomly. 
We present representative results for N — 10 for Landau- 
Stuart and Rossler oscillators in middle and right pan- 
els respectively in Fig. [5] The rotations of individual 
oscillators are indicated in Fig. [5Ja,b,c) where open 
and filled circles correspond to the clockwise and anti- 
clockwise rotations respectively in these particular real- 
izations. Three different topologies are taken: (i) near- 
est neighbor coupling with periodic boundary conditions 



(K = 2)~Fig. EJa); (ii) global coupling where each oscil- 
lator is coupled with all others (K = N — 1) - see Fig. 
[SJb); (iii) each oscillator is coupled to three (K = 3) 
random nodes — see Fig. [5jc) which realizes a small 
world topology [la ]. The coupling functions H for these 
network of oscillators are the same as in Figs. [D and 0] 
for Landau- Stuart and Rossler oscillators respectively. In 
Figs. [SJd) to (i) the largest few Lyapunov exponents are 
plotted as a function of the coupling strength in the three 
cases. The inset figures in Fig. [5Jd-i) show the relative 
phase xa vs xb (black curves) and yAVSys (red curves) 
for counter-rotating oscillators between the marked os- 
cillators A and B in (a) verifying the existence of mixed- 
synchronization for both systems. The occurrence of AD 
is also observed in networks of Landau-Stuart oscillators, 
as shown by negative Lyapunov exponents in the middle 
panel. 

In summary, we have observed the universal occurrence 
of mixed-synchronization in counter-rotating oscillators: 
some variables are synchronized in-phase, while others 
are out-of-phase. The phenomenon of amplitude death is 
also observed in some cases, as for example, the Landau- 
Stuart limit cycle oscillators studied here. These results 
are also found to persist in networks of oscillators of dif- 
ferent topologies. We have also checked that these results 
persist under parameter mismatch |l4j, and believe that 
they are generally applicable and should be observable 
in experiments. Thus the phenomena observed here will 
occur in other coupled systems as well. 

In this paper we have considered only linear diffusive 
coupling, and other types of interactions such as one-way, 
conjugate [ll[, an d nonlinear couplings 12 1 are still need 
to be explored. Dynamical phenomena such as riddling, 
hysteresis, phase locking, phase-shifting, phase-flip [l|-|6| 
etc. which have been seen in co-rotating coupled systems 
should also be explored in counter-rotating coupled os- 
cillators. 
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